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ABSTRACT 

An  experimental  9th  grade  computer  science  syllabus 
is  proposed.  The  syllabus  would  include  the  technical  information 
needed  for  controlling  and  programing  the  computer  in  a number  of 
modes  and  would  preview  some  of  the  areas  covered  in  the  high  school 
curriculum.  A sample  module  of  a topic  not  normally  taught  in  high 
school — distance  and  error- correcting  codes — is  presented.  Two 
student-authored  programs  also  illustrate  the  ways  in  which  the 
computer  can  explore  areas  outside  the  normal  curriculum — a plotter 
routine  for  making  artistic  patterns  and  a program  to  use  the  Monte 
Carlo  method  to  calculate  the  area  under  a curve.  (JY) 
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Computer  Science  Modules 

Our  recent  experience  indicates  that  there  is  a lot  of  potential  for 
getting  youngsters  in  the  9^,  8^,  (or  sometimes  7^)  grade  excited  about 
exploring  the  world  of  learning  with  computers.  The  potential  is  magnified 
several  times  when  the  student  controls  the  full  power  of  the  computer;  when 
it  is  used  by  him,  not  on  him.  This  means  that  we  must  provide  the  oppor- 
tunity for  these  students  to  master  control  of  the  machine  at  an  early  a.ge. 

The  implications  for  increased  sophistication  of  use  in  the  upper  grades  is 
obvious. 

To  this  end  we  are  proposing  that  an  experimental  9**1  grade  "Computer 
Science"  syllabus  be  prepared  for  testing  next  fall.  Such  a syllabus  would 
include  the  technical  information  needed  for  controlling  and  programming  the 
computer  in  a number  of  modes,  with  applications  that  "previewed"  some  of 
the  things  the  student  could  learn  more  about  in  his  high  school  years  ahead. 
For  many  students  this  could  mean  a new  sense  of  purpose  in  working  at 
their  education. 

The  computer  science  syllabus  would  also  include  topics  not  normally 
taught  in  high  school  which  have  a "fun11  characteristic  to  them.  A sample 
module  illustrating  such  a topic  (DISTANCE  AND  ERROR- CORRECTING 
CODES)  is  enclosed. 


Two  Student- Authored  Programs 

The  two  programs  enclosed  also  illustrate  topics  outside  the  usual  cur 
riculum.  The  first  is  a plotter  routine  for  making  artistic  patterns  written 
by  Chris  Van  Sickle,  an  8^  grade  student  who  "sneaked"  into  our  program. 
The  second  (Monte  Carlo  Integration)  is  by  Mike  Kaufmann  who  is  a senior 
at  Allderdice. 
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Information  Retrieval 


A preliminary  list  of  module  areas  that  come  under  the  heading  of  in- 
formation retrieval  is  included  with  this  newsletter.  Some  of  these  modules 
could  supplement  mathematics  courses.  They  all  can  be  classified  under 
computer  science.  Students  who  master  these  techniques  can  apply  them  to 
any  subject. 

^Supported  in  part  by  NSF  grant  GJ  1077 

1 


( 


2 


PRELIMINARY  LIST  OF  MODULE  AREAS  UTILIZING 
INFORMATION  RETRIEVAL 


1.  Set  Theory:  Basic  Concepts 

a.  Definition 

b.  Examples 

c.  Modification  of  Set  (updating) 

2.  Operations  on  Sets 

a.  Union 
Intersection 
Complement 

Commutative,  associative  and  distributive  laws 

b.  Logical  interpretation  of  operations 

3.  Information  Retrieval  Techniques 

a.  Description 

b.  Examples 

c.  Building  and  updating  file 

(1)  selection  of  material 

(2)  representation 

(3)  structuring  ^ 

(4)  set  interpretation 

d.  Searching 

(1)  constructing  query 

(2)  query  analysis 

(3)  logical  operations 

(4)  iterating  the  search 

4.  Information  Retrieval:  Advanced  Concepts 

a.  Relevance 

b.  False  drop  analysis 

c.  Modifying  representation 

d.  Syntactic  and  semantic  problems 

e.  Fact,  reference,  and  information,  retrieval 

5.  Guided  Application  of  Information  Retrieval 

a.  Select  subject  matter  (e.g.  Apollo  Moon  Program) 

b.  Student  builds  I.R.  system  using  techniques  given  in 
3 above . 

c.  Student  uses  retrieval  system  to  study  some  pertinent 
problem  in  the  subject  under  consideration 

6.  New  I.R.  systems 

Student  is  to  write  and  demonstrate  his  own  fact  retriev- 
al system.  (Some  subjects  might  be  professional  football 
scores,  prices  of  a given  stock  over  the  past  30  years, 
data  on  the  planets,  super-market  prices,  areas  suggested 
by  the  school  librarian,  etc.) 


2 


• v I 


Distance  and  Error- 
Correcting  Codes 
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A binary  code  of  length  N is  a string  of  N 0's  or  l's. 
For  example,  if  N = 3,  all  the  possible  binary  codes  are  000, 
001,  010,  011,  100,  101,  110,  and  111.  We  speak  of  these  as 
3 BIT  codes  (1  and  0 are  called  BITS) . 

These  codes  could  be  used  to  represent  eight  objects  of  any 
sort  the  members  of  a musical  octet,  the  digits  0,1, 2, 3, 4, 5, 6, 7 
in  a computer,  or  the  letters  A,B,C,D,E,F,G,H. 


Now  for  some  intrigue 

Let's  assume  that  we  wish  to  assign  binary  code  names  to 
the  agents  of  STICK  (Society  to  Increase  Contact  for  Keeps) , an 
international  rinc  of  glue  thieves.  Suppose  we  only  have  two 
agents  but  eight  codes.  Question:  Can  we  assign  codes  so  that: 

(a)  The  computer  will  check  code  authenticity  without 
knowing  the  correct  codes. 

(b)  The  computer  can  give  the  correct  code  even  though 
the  agent  has  deliberately  changed  one  BIT  (to  throw 
off  eavesdroppers) . 

To  see  how  codes  can  be  assigned  to  make  this  possible, 
let's  place  the  codes  at  the  vertices  of  a cube. 


H . ^°  k®  more  precise,  we  should  call  the  above  figure  a 
3-dimensional  cube".  A picture  of  a "4-dimensional 'cube" 

(which  has  2 = 16  vertices)  is  shown  on  the  cover.  Thus  we  can 
associate  a unique  four-bit  code  with  each  vertex  of  a 4-D  cube. 
Can  you  generalize  this  statement? 


4 


We  will  define  the  number  of  BITS  by  which  two  codes  differ 
as  the  "DISTANCE" (=d)  between  these  codes.  Thus,  for  example, 
the  distance  between  001  and  101  is  d = 1,  the  distance  between 
001  and  111  is  d = 2,  and  tiie  distance  between  001  and  110  is 
d = 3.  Math  students:  Is  this  a legal  use  of  the  word  distance? 
Notice  that  our  picture  has  been  drawn  so  that  "distance"  between 
codes  corresponds  to  the  number  of  edges  of  the  cube  you  would  Lave 
to  walk  along  to  get  from  one  vertex  to  the  other. 

Let's  assign  our  two  authentic  agents  the  codes  001  and  110 
(which  are  a distance  of  three  from  each  other) . Now’  suppose  one 
agent  walks  up  to  another  and  says  my  code  is  101. 

(a)  How  can  we  tell  if  it  is  an  authentic  code?  One  way 
would  be  to  simply  compare  it  to  the  list  of  authentic 
cod€?si  However,  if  there  were  to  be  very  many 
codes,  such  a search  of  the  authentic  list  would  be 
time  consuming.  Besides,  we  don't  want  this  authentic 
list  stored  in  too  many  places!  There  is  another  way 
to  check  authenticity. 

In  our  example  the  two  authentic  codes  have  the  pro- 
perty that  if  we  add  the  first  and  third  BITS  of  the 
code  we  get  1,  and  this  is  also  true  if  we  add  the 
second  and  third  BITS. 


110 
A 


sum 
sum=l 


* 


For  all  other  codes  this  is  false. 

NOTE:  0+1=1  and  1 + 0 = 1 in  binary  arithmetic. 

ALSO  NOTE:  0+0=0,  but  1+1=0  (with  carry  of  1) 

FURTHER:  0*1=0  1*0=0  0*0=0  and  1*1=1 

Thus  the  code  1 0 1 our  agent  gave  is  not  authentic. 

Y\ 

Since:  V-l  + 1 = 0^ 
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(b)  Suppose  our  agent  deliberately  changes  one  BIT  in  his 
code  when  giving  it  verbally.  Applying  the  above  rule 
will  detect  the  error,  but  can  we  figure  out  what  the 
correct  code  should  have  been? 

We  can  see  the  answer  from  our  diagram.  An  authentic 
code  which  has  only  one  BIT  changed  is  distance  d = 1 
from  the  original  correct  code,  but  distance  d = 2 from 
the  other  correct  code.  Thus  101  has  to  be  corrected 
back  to  001,  not  to  110, 

Try  to  develop  an  algorithm  for  making  corrections  in 
codes  which  have  1 BIT  in  error : Here  is  how  you  might 
reason: 

Let's  call  the  3 BITS:  Bl,  B2,  and  B3. 

FOR  101  Bl  + B3  = 0 WRONG 

B2  + B3  = 1 RIGHT 

.*.  Change  Bl,  Correct  code  is:  001 

FOR  111  Bl  + B3  = 0 WRONG 

B2  + B3  = 0 WRONG 

.*.  Change  B3,  Correct  code  is:  110  etc. 


Problems:  Write  programs  for  your  wrist  computer  to  handle 

the  following: 

1 


1. 


2. 

3. 


4. 


The  Two-Agent  Problem. 

Authentic  Codes:  101  and  010 


INPUT:  Any  3 BIT  code  which  is  either  an  authentic 

code,  or  which  contains  an  error  in  1 BIT. 


OUTPUT : 


The  Message: 

OR: 


"AUTHENTIC  CODE" 
"CODE  IN  ERROR 
CODE  SHOULD  BE 


Can  use  of  a 4 BIT  code  (see  picture  on  cover)  permit  add- 
itional outputs  for  the  above  "Two-Agent"  analysis  program? 


Four-Agent  Problem. 
Authentic  Codes  are: 


INPUT : Any  code 

OUTPUT:  The  Message: 

OR: 

OR: 


SMITH  00000 

BOND  11100 

SPIRO  00111 

JONES  11011 


"AUTHENTIC  CODE" 

"1  BIT  ERROR — CORRECT  CODE  IS 
"ERROR  >=2  BITS— DOUBLE  AGENT11’ 


Here  is  a set  of  six  BIT  codes  to  play  with:  000000,  000111, 
111000,  110110,  011011,  101101.  (NOTE:  d>=3  for  any  two  of 
these  codes.) 
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Sample  Solution — Problem  1 

>L I STNH 

10  SC  1 > - 0 
£0  SC2>  = 1 

8 5 PR* "TWO- AGENT  PROBLEM" 

9 0 "TO  END  PROGRAM*  TYPE  THE  CHARACTERS  END  WHEN  ASKED  FOR  A CODE 

• • 

9 5 PR."  " 

100  PR.  "ENTER  A 3-BIT  CODE." 

105  INPUT  B$ 

106  IF  B 35  = "END"  GOTO  185 
1 10  GO SUB  200 

115  SC  3)  = SI 
120  SC  4)  = S2 

125  IF  CSC1)  = SC3)  AND  SC2)=SC4>>  GOTO  IPO 

126  IF  CSC1)  it  SC  3)  AND  SC2)0SC4>>  GOTO  140 

129  X = 1 -SC  3) 

130  IF  CSC1)  = X)  GOTO  150 

1 34  X=  1 - SC  4) 

135  IF  C SC  2)  = X)  GOTO  160 

140  B3S  = RIGHTC  STRC  j -VALCB335) ) il) 

1 45  GOTO  165 

150  B 1 = RIGHTC  STRC 1 - VAL  C B 1 $) ) * 1) 

155  GOTO  165 

160  B2£  = RIGHTC  STRC 1 - VALCB2S))*  1) 

165  PRINT"  CODE  IN  ERROR" 

170  PRINT"  CODE  SHOULD  BE  "S  B 1 S+B2S  + B3S 
1 72  GOTO  95 
1 75  STOP 

180  PRINT"  AUTHENTIC  CODE" 

182  GOTO  95 
185  STOP 
190  END 

200  BIS  = LEFTCBS, 1 ) 

205  X$  = RIGHTCB$,2> 

2 1C  B2S  = LEFTCXS*  1 ) 

2 15  B3S  = RIGHTCBS* 1 ) 

220  SI  = VALC  BIS)  + VALCB3S) 

225  S2  = VALCB2S)  + VALCB3S) 

230  IF  SI  = 2 LET  SI  = 0 
235  IF  S2  = 2 LET  S2  = 0 
240  RETURN 
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Problem  1 (cont. ) 


>RUN 

TW0 -AGENT  PROBLEM 

T®  END  PR0GRAM,  TYPE  THE  CHARACTERS  END  WHEN  ASKED  F0R  A C0DE 

ENTER  A 3-BIT  C0DE. 

?000 

C0DE  IN  ERR0R 
C0DE  SH0ULD  BE  010 

ENTER  A 3-BIT  CODE. 

? 001 

C0DE  IN  ERROR 
CODE  SHOULD  BE  101 

ENTER  A 3-BIT  CODE. 

?010 

AUTHENTIC  CODE 

ENTER  A 3-BIT  CODE' 

? 01  1 

C0DE  IN  ERROR 
CODE  SHOULD  PE  010 

ENTER  A 3-BIT  CODE. 

? 100 

CODE  IN  ERROR 
CODE  SHOULD  BE  101 

ENTER  A 3-BIT  CODE. 

? 101 

AUTHENTIC  CODE 

ENTER  A 3-BIT  CODE. 

? 1 10 

CODE  IN  ERROR 
CODE  SHOULD  BE  010 

ENTER  A 3-BIT  CODE. 

? 1 1 1 

CODE  IN  ERROR 
CODE  SHOULD  BE  101 


ENTER  A 3-BIT  CODE 
? END 
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Sample  Solution — Problem  3 
>LI STNH 


PR 

VAR  = ZER0 
PR • "ENTER  A 
INFUT  B$ 

IF  B$  = "END*  GOTO 
BIS  = LEFTCBS* 1) 

XS  = RI GHTC  B$*  4) 
B2S  - LEFTC XS* 1 ) 

XS  = RI GHTC  B$*  3) 
B3S  “ LEFTC  XS*  1 ) 

X$  = RIGHTCB$*2) 
B4$  = LEFTC  XS* 1 ) 
B5S  = RI GHTC  X $*  1 ) 


5-BIT  BINARY  CODE" 


250 


SC  1 ) 

SC  2 ) 

SC  3 ) 

SC  4) 

F0R  I - 
IF  SCI) 
NEXT  I 
IF  SCI) 


VALC  BIS) 
VALCB4S) 
VAL  C B 1 S ) 
VALCB2S) 
1 T0  4 
= 2 LET 


+ VALCB2S) 
+ VALCB5S) 
+ VALCB^'S) 
+ VALCB5S) 

SCI)  = 0 


IF  SC  2) 

IF  FI  = 

F2  - 1 
IF  VALCB3S) 
IF  VALCB3S) 
K = FI  + F2 


ft  0 LET  FI  = 
= 0 G0T0  125 
1 G0T0  240 


5 PR.  "FOUR-AGENT  PROBLEM" 
10  PR. "T0  END  PROGRAM*  TYPE 

| p DD  _ »•  »* 

1 3 
1 5 
20 
25 
30 
35 
40 
45 
50 
55 

6 0 
65 
70 
7 5 
80 
85 
90 
95 
1 00 
1 05 
1 10 
1 15 
1 20 
1 25 
1 30 
1 35 
1 40 
1 45 
1 50 
1 55 
1 60 
1 65 
1 70 

1 75 
1 80 
185 
190 
195 
200 
205 

2 10 
215 
220 
225 
2 30 
235 
2 40 
245 
250 
260 


THE  CHARACTERS  END  WHEN  ASKED  F0R  A CODE* 


ft  SC  3)  LET  F3  = 1 
ft  SC  4)  LET  F 4 = 1 
+ F3  +.  F 4 
G0T0  230 
G0T0  24C 

F3 ) ft  2 G0T0  165 

STRC1  -VALCB1S) )* 1 ) 


ft  2 G0T0  180 

STRC1  - VALC  B2S ) ) * 1 > 


IF  K = 0 
IF  K ft  2 
IF  CF1  + 

BIS  = RIGHTC 
G0T0  215 
IF  CF1  + F4) 

B2S  = RIGHTC 
GOTO  215 
IF  CF3  + F4) 

B3S  = RIGHTC 
G0T0  215 
IF  CF2  + F3 ) 

B4S  = RIGHTC 
G0T0  215 
B5S  =RI GHTC 
PR."  " 

1 BIT  ERROR CORRECT 

12 

AUTHENTIC 
12 

ERROR  = 2 

12 


it  2 GOTO  195 
STRC1  - VAL C B3S ) ) * 1 ) 

#2  G0T0  210 
STRC 1 - VALC  B4S ) ) * 1 ) 

STRC 1 - VALCB5S) )* 1 ) 


PR." 

G0T0 

PR." 

G0T0 

PR.%: 

GOTO 

STOP 

END 


CODE  IS  "SB! S+B2S+B3S+B4S+B5S 


CODE’ 


BITS  DOUBLE  AGENT!" 


O 

ERIC 


9 


Problem  3 (cont.) 
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>RUN 

FOUR-AGENT  PROBLEM 

T0  END  PROGRAM*  TYPE  THE  CHARACTERS  END  WHEN  ASKED  F0R  A C0DE 

ENTER  A 5-BIT  BINARY  C0DE 
? 00000 

AUTHENTIC  C0DE 

ENTER  A 5-BIT  BINARY  C0DE 
7 10000 

1 BIT  ERR0R CORRECT  CODE  IS  00000  - 

ENTER  A 5-BIT  BINARY  CODE 
700100 

1 BIT  ERROR CORRECT  CODE  IS  00000 

ENTER  A 5-BIT  BINARY  CODE 
700010 

1 BIT  ERROR CORRECT  CODE  IS  00000 

ENTER  A 5-BIT  BINARY  CODE 
711011 

AUTHENTIC  CODE 

ENTER  A 5-BIT  BINARY  CODE 
7 11010 

1 BIT  ERROR CORRECT  CODE  IS  11011 

ENTER  A 5-BIT  BINARY  CODE 
7 10011 

1 BIT  ERROR CORRECT  CODE  IS  11011 

ENTER  A 5-BIT  BINARY  CODE 
710101 

ERROR  = 2 BITS  DOUBLE  AGENT! 

ENTER  A 5-BIT  BINARY  CODE 
710110 

ERROR  = 2 BITS  DOUBLE  AGENT! 

ENTER  A 5-BIT  BINARY  CODE 
7 END 
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THE  MONTE  CARLO  METHOD 


—A  FUN  STATISTICAL  APPROACH  TO  CALCULATING  THE  AREA  OF  , 
GEOMETRICAL  FIGURES  OR  THE  AREA  UNDER  A CURVE 

(by  Michael  Kaufman,  Grade  12) 

Method ; 

Set  off  a rectangle  which  includes  the  area  to  be  calculated. 
Using  the  candom  number  generator,  try  out  many  points  included 
in  the  rectangle,  to  see  if  they  are  also  included  in  the  figure. 
The  number  of  points  included  in  the  figure,  divided  by  the 
total  number  of  points  tried,  times  the  area  of  the  rectangle, 
will  approximate  the  unknown  area. 


Problem: 


Find  the  area  of  a quarter  of  a unit  circle  using  the  Monte 


(0,1 VA 


Carlo  method. 

Set  off  rectangle  (square)  ABCD 
so  that  it  includes  the  quarter  cir- 
cle whose  area  is  to  be  calculated. 

We  will  'shoot'  500  points  at  the 
square.  The  points  each  have  co- 
ordinates (x,y)  where  both  x and 
y are  between  0 and  1 (this  puts 
the  points  within  the  square  ABCD) . 

The  (x,y)  coordinates  will  be  found 
using  the  random  number  generator. 

The  equation  of  a unit  circle  is 

2 2 2 2 t 

fore  x + y = 1 , y = 1 - x,  y=  /1-x* 2.  (The  /l  - x2  eliminates 
all  negative  values  of  y so  that  this  equation  describes  a semi- 
circle. Since  we  are  taking  x between  0 and  1,  we  will  be 
finding  one-half  of  that  semicircle.  If  we  had  taken  -1  < x < 1, 
we  would  have  a semicircle.) 

The  number  of  points  (x,y)  such  that 


x2+  V2=  r2. 


r = 1 . There- 


< A - x: 


divided 


by  the  total  number  of  points,  times  the  area  of  the  square, 
equals  the  area  of  the  quarter  circle  (it  should  be  ^ 


7T  * 2 

* r 


r = 1,  therefore  AREA  = - = .785). 
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Let's  see  what  the  computer  gets: 

5 T«-A«-Y«-N«-M«-0  (or— 5 VAR=ZERO) 

10  LET  N=500 

20  FOR  1=1  TO  N 
30  LET  X=NUM(N)/N 

40  LET  Y=NUM(N)/N 

50  IF  Y<SQRT (l-X+2 ) LET  T=T+1 

60  LET  M=M+1 

70  NEXT  I 
80  LET  A=  (T/M)  *1 

90  PR.  "FOR  " : I : " POINTS,  AREA= " : A 
100  END 

>RUN 

FOR  500  POINTS,  AREA=.756 


Multiple  assignment  setting 
T,  A,  Y,  N,  M all  = 0. 

Number  of  points  we're  trying. 

Generates  500  possible  val- 
ues of  x for  0 < x < 1. 

Same  for  y. 

If  point  is  within  quarter  cir- 
cle, increase  counter  by  one. 

Counts  total  number  of  points 

Area  = points  within  the  quar- 
ter circle/total  points,  times 
the  area  of  the  square. 


The  more  points  (up  to  a limit) 
the  more  accurate  the  area. 


Caution : 

Pay  particular  attention  to  the  range  of  x and  y.  For 
instance,  can  you  see  why  steps  30  and  40  in  the  above  program 
would  be  meaningless  for  a circle  of  radius  2?  The  computer 
would  respond:  FOR  500  POINTS,  AREA= 1.000  Why? 

Because  the  rectangle  is  totally  within  the  circle,  so  the 
computer  is  really  calculating  the  area  of  a square  of  side  1 
which  happens  to  be  = 1.000.  For  a circle  of  radius  2,  the  range 
for  x and  y must  be  [0,2]. 

Hint : 

A useful  formula  which  helps  to  set  the  range  on  the  random 
numbers  generated  is:  FOR  A<B , (B-A) *R+A  , where  A and  B are 

the  needed  limits  and  R is  a random  number  generator,  generating 
values  between  0 and  1.  For  instance,  we  want  the  range  for  a 
certain  calculation  to  be  7 and  11.  The  formula  we  use  is:  x = 
(11-7) *NUM(N)/N+7.  For  R = 0,  x = 7;  for  R = 1,  x = 11. 

' v~ 

R 

Your  Turn: 

1.  Use  1,000  points  to  obtain  a better  approximation  of  the  area 
of  the  quarter  qircle. 


12 


MONTE  CARLO 
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2.  Modify  the  program  on  page  2 so  that  the  user  can  input  any 
radius.  Don't  forget  the  limits  of  x and  yl  Also,  Have 
the  program  give  the  area  of  the  entire  circle. 

3.  Find  the  area  under  the  curve  y = x for  x between 
0 and  10.  Use  an  appropriate  set  of  limits  for  y. 

4.  Find  the  area  under  the  curve  y = sin  (x)  for  x between 
0 and  it. 

3 2 

5.  Find  the  area  under  the  curve  y = x + 2x  +3  for  0 < x < 1. 
Use  x between  0 and  1,  and  y between  0 and  6.  Then  ask 
your  teacher  to  use  integral  calculus  to  get  an  exact  area 
under  the  same  curve.  Compare  the  two  values. 

6.  Now  show  your  teacher  that  you  can  easily  find  the  area  under 
the  curve:  y = x2/ (x2+  32 ) 3 , for  0 < x < 1,  while  he  may 
find  it  not  quite  so  easily.  Hint  to  the  student:  Use  x 
between  0 and  1,  and  y between  0 and  31.62278. 

Hint  to  the  teacher:  Use  0 and  1 as  the  limits  of  integration. 

/ x2/.(x2+  a2  ) 3dx  - — / (x2+  a2  ) 5 - — ~/(x2  + a2)3 

J 6 24 

4 6 

- ^p~Vx2+  a2  - Yg’log  (x2+ /x2+  a2) 

Good  luck  I 

7.  Using  a given  curve,  plot  your  accuracy  versus  the  number  of 
points  used  in  the  computer  run.  For  instance,  use  the  circle 
program  on  page  2.  Try  10  points,  100  points,  200  points, 

300  points,  etc.,  and  plot  on  a graph  the  increasing  accuracy 
versus  the  number  of  points  per  trial. 

8.  See  if  you  can  think  of  some  more  interesting  applications 
of  the  Monte  Carlo  method. 

Optional : 

1.  Choose  a very  small  area  (such  as  .0001  by  .0001)  and  see  if 
the  Monte  Carlo  method  is  any  better  for  a small  area. 

2.  Graph  increasing  area  versus  accuracy.  (Note:  accuracy 
should  be  calculated  as  the  percentage  difference  between 
the  Monte  Carlo  answer  and  the  "correct"  answer.) 

3.  Use  the  H-P  Plotter  to  graph  your  curve,  and  then  to  plot  each 
of  the  random  points.  Don't  forget  to  also  plot  the  rec- 
tangular region.  You  will  get  a very  good  idea  of  how  Monte 
Carlo  works  by  watching  this  process. 

o 

ERIC 
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THE  MONTE  CARLO  METHOD:  TEACHER'S  GUIDE 


The  main  problem  that  students  might  have  with  this  method 
is  understanding  how  the  ranges  on  x and  y are  set.  What 
should  be  emphasized  is  that  the  maximum  value  of  the  ranges  are 
the  dimensions  of  the  rectangle  placed  around  the  figure. 

Another  weak  point  is  that  students  might  not  understand  why 
the  fraction  T/M  (the  number  of  points  in  the  figure/the  total 
number  of  points  tried)  should  be  multiplied  by  the  area  of  the 
rectangle.  The  fraction  T/M  represents  the  fraction  of  the 
rectangle  taken  up  by  the  figure.  This  fraction  of  the  rectangle 
multiplied  by  the  area  of  the  rectangle  will  give  the  area  of  the 
figure. 


Students  could  be  introduced  to  the  integral  calculus  using 
this  module.  They  should  see  that  where  the  Monte  Carlo  method 
provides  a very  approximate  estimate  of  the  area  under  a curve, 
integration  is  defined  to  be  the  exact  area  under  the  curve  be- 
cause it  takes  the  limit  of  an  infinite  number  of  dx's;  similarly 
we  could  take  an  infinite  number  of  points  and  come  out  with  a 
precise  value. 


Answers  to  'Your  Turn': 


1. 

Modify  step  10  to  read: 

LET  N=1000 

2. 

3. 

Modify  program  on  page  2: 
Ranges : 0 < y < 10 , 0 < x < 10 

8 

30 

40 

30 

LET  X=10*NUM(N)/N 

50 

40 

LET  Y=10*NUM(N)/N 

80 

50 

IF  Y<X  LET  T=T+i 

95 

80 

LET  A=T/M*100 

AREA=50 

10 

4. 

Modify  program  on  page  2: 

Ranges:  0 < y < 1 , 0 < x < 3 . 14159265 

30 

LET  X=  3.14159265 *NUM (N ) /N 

rr 

50 

IF  Y<SIN  (X)  LET  T=T+1 

0 ■ 

80 

LET  A=T/M*3. 14159265 

A 

AREA=2 . 0 

40 

50 

6. 

Modify  program  on  page  2: 

80 

Modify  program  on  page  2: 

Ranges : 0<y<R,  0<x<R. 

PR.  "INPUT  RADIUS" : INPUT  R 
LET  X=R*NUM(N)/N 
LET  Y=R*NUM(N)/N 
IF  Y SQRT (R  2-X  2)  LET  T=T+1 
LET  A=T/M*R  2 

PR.  "AREA  OF  ENTIRE  CIRCLE=" : 4*A 

Might  also  include: 

PR.  "INPUT  N": INPUT  N 


Modify  program  on  page  2: 
Ranges : 0 < y < 6 , 0 < x < 1 

LET  Y~6*NUM(N)/N 

IF  Y< (Xi3+2*Xi2+3)  LET  T=T+1 

LET  A=T/M*  6 . 0 

AREA=3.91666  (by  integration 

and  slide  rule) 


40  LET  Y=31. 62278*NUM(N)/N 
(31 . 62278=SORT (1000)  ) 


50  IF  Y< ( (Xi2) *SQRT ( (Xi2+9) i3) ) LET  T=T+1 
80  LET  A=T/M*31. 62278 
ARE A=  ?????????? 


Note:  Use  a FOR-NEXT  loop  to  automatically  increase  N in  problem  7. 


□ram  — innLQzn] 
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7 OPEN  //  FOR  0UTPUT  AS  3 
l 0 K = 0 

12  B=NUMC100> 

1 4 IF  B<25  GOTO  12 
1 6 X=NUM(9999> 

18  Y=MUM(8000) 

19  IF  Y-sJOOO  GOTO  18 

20  Z=NUM< 9999 ) 

22  M=INT<CX+Z>/2> 

24  N=INTC ABS<X-Z>/2> 

26  IF  ARS<  X-Z ) < 1 000  0R  ABSCX-Z ) > 5000  G0T0  16 
30  PRINT  0N  3 1 "PLTP" 

32  PRINT  0N  31  XJY 

34  PRINT  0N  3J  "PLTL" 

36  PRINT  0N  3:  Z;Y 

38  PRINT  0N  31  "PLTP" 

40  F0R  I =B  BY  B 

41  IF  I >M  0R  Z~I<0  0R  X-I<0  0R  Y-I<0  0R  Z+I>9999  0R  X+I>9999  0R  Y+I>99 
99  G0T0  62 

42  IF  X<Z#  X =X  + 1 ELSE  X = X-I 

43  IF  X<Z*  Z=Z-I  ELSE  Z=Z+I 

44  PRINT  0N  3:  X;Y  + I 

45  PRINT  0N  3S  "PLTL" 

46  PRINT  0N  3s  ZiY-I 

47  PRINT  0N  3!  "PLTP" 

48  PRINT  0N  3:  ZlY+I 

49  PRINT  0N  3:  "PLTL" 

50  PRINT  0N  3:  XlY-I 

5 1 PRINT  0N  3:  "PLTP” 

52  NEXT  I 

62  L=L* 1 * CIF  L>3*  G0T0  663 

63  TYPE  "ONE  D0NE" 

64  G0T0  12 

66  PRINT  0N  3S  "PLTT" 

67  K=K+ 1 * L = 0#  CIF  K>3#  G0T0  763 
7 4 G0T0  12 
7 6 CL0SE  3 
7 8 END 
> RUN 

0 PEN I NG  FILE  3,  SYMBOLIC  OUTPUT » NAMES  /STARS/ 

NEW  FILE? 

ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 
ONE  DONE 


Listing  of  a plotter  program 
written  by  Chris  Van  Sickle 
(Grade  8) . Part  of  the 
pattern  produced  on  the 
plotter  is  reproduced  on 
the  cover . 


> -COPY 

FILE.  NAMES/STARS/ 
TO  s TPT 
PLTP 

7047  3483 

etc-  . . 
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